3 Y NMR Probe of Zn Induced Local Magnetism in YBa 2 (Cu! _ y Zn y ) 3 
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We present detailed data and analysis of the effects of Zn substitution on the planar Cu site in 
YBa2Cu306+a; (YBCOg+i) as evidenced from our 89 Y NMR measurements on oriented powders. 
For x « 1 we find additional NMR lines which are associated with the Zn substitution. From our 
data on the intensities and temperature dependence of the shift, width, and spin-lattice relaxation 
rate of these resonances, we conclude that the spinless Zn 3d 10 state induces local moments on the 
near-neighbour (nn) Cu atoms. Additionally, we conjecture that the local moments actually extend 
to the farther Cu atoms with the magnetization alternating in sign at subsequent nn sites. We 
show that this analysis is compatible with ESR data taken on dilute Gd doped (on the Y site) and 
on neutron scattering data reported recently on Zn substituted YBC06+a;. For optimally doped 
compounds 89 Y nn resonances are not detected, but a large T-dependent contribution to the 89 Y 
NMR linewidth is evidenced and is also attributed to the occurence of a weak induced local moment 
near the Zn. These results are compatible with macroscopic magnetic measurements performed on 
YBCOs+a; samples prepared specifically in order to minimize the content of impurity phases. We 
find significant differences between the present results on the underdoped YBCO6+1 samples and 
27 Al NMR data taken on Al 3+ substituted on the Cu site in optimally doped La2Cu04. Further 
experimental work is needed to clarify the detailed evolution of the impurity induced magnetism 
with hole content in the cuprates. 
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I. INTRODUCTION 

It is now experimentally well established that the 
Cu02 planes are responsible for the magnetic and su- 
perconducting properties of the cuprates. However the 
interconnection between these two properties is still an 
essential but unanswered question. Understanding the 
normal state of the cuprates is still a prerequisite for any 
theoretical approach to the microscopic origin of High 
Temperature Superconductivity. In the recent past, con- 
siderable interest has been aroused due to the detection 
of a pseudo-gap in the spin excitation spectrum of the 
cuprates for underdoped materials (the word pseudo is 
prefixed because although a strong decrease in the inten- 
sity of excited states is detected well above the super- 
conducting transition temperature T c , a real gap is only 
detected below T c ). The first indications of a pseudo-gap 
were provided in the microscopic NMR measurements of 
the susceptibility of the Cu0 2 planes. 89 Y NMR shifts 
in YBCOg+o; of Alloul et al. Q| were found to decrease 
markedly with decreasing T. The large decrease of the 
static susceptibility was interpreted to be due to an open- 
ing of a pseudo-gap in the homogeneous q = excitations 
of the system. In the underdoped (x < 1 for YBCOe+z) 
high-T c cuprates, a similar decrease of the spin-lattice 
relaxation rate 63 Cu \jT\T O, which is dominated by 
the imaginary part of the susceptibility at the AF wave 



vector q = (jr, w) was also indicative of a pseudo-gap in 
the spin excitations at this wave vector Q . 

The inelastic neutron scattering experiments which fol- 
lowed , clearly confirmed the existence of this pseudo 
spin-gap at (71", n), and allowed measurements of the fre- 
quency dependence of the excitations. The temperatures 
at which these two pseudo-gaps begin to open are found 
to be different, and it is not clear at present whether 
they signal different cross-overs between distinct states 
or a single cross-over phenomenon. The latter case 
would imply a wave-vector dependence of the pseudo- 
gap. Presently, the existence of a pseudogap for the un- 
derdoped high-T c superconductors has been detected by 
many techniques such as transport, photoemission, etc. 
While various explanations are proposed for the pseudo- 
gaps, it is believed that the essential physics of the nor- 
mal state (and perhaps the superconducting state), of 
the cuprates might be linked to it. 

In order to better characterise the properties of the 
cuprates it has appeared quite important to understand 
their modifications due to impurities or disorder. Atomic 
substitutions on the planar Cu site are naturally found 
to be the most detrimental to superconductivity, while 
modification in the charge reservoir chains mainly yield 
changes in hole doping. For such studies the YBCOg+z 
system is particularly suitable, as the variation of impu- 
rity induced magnetism with hole doping can be studied 
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by merely changing the oxygen content x for a given im- 
purity content. In classical superconductors, T c is mainly 
affected by magnetic impurity substitutions. In cuprates 
it has been shown that even a non-magnetic impurity like 
Zn (3d 10 ), which substitutes on the Cu site of the CuC>2 
plane strongly decreases the superconducting transition 
temperature T c (about 10.6 K/% Zn for x — 1). It has 
also been anticipated || and then shown experimentally 
that although Zn itself is non-magnetic, it induces a mod- 
ification of the magnetic properties of the correlated spin 
system of the CuC>2 planes |7j . Using 89 Y NMR we have 
further shown, in the preliminary report of the present 
work ||, that local magnetic moments are induced on 
the nn Cu of the Zn substituent in the CuC>2 plane. Two 
important result have been demonstrated: 

i) the q — pseudo-gap was found unaffected by Zn 
even when T c is reduced to zero for YBCOg.6, 

ii) The magnitude of the induced local moment is 
strongly dependent on the carrier concentration || . 

Since our reports, other experimental evidence by 

NMR in YBCO 0,0, in YBa 2 Cu 4 8 (1248) ||, in 
La 2 Cu0 4 @, or ESR in YBCO and 1248 Q, have con- 
firmed that the occurrence of local moment induced by 
non-magnetic impurities on the Cu sites is a general prop- 
erty of cuprates. The local moments have been observed 
as well in macroscopic bulk susceptibility data |l7| . 
The Zn-induced modifications of the magnetic excitations 
both in the superconducting and the normal state have 
been studied by neutro n scattering 18 19(|. Also, 
trical transport 



21|, and thermal properties [[22j of 
substituted high-T c cuprates have been investigated. 

However some studies have concluded that the suscep- 
tibility near the Zn does not exhibit a Curie behaviour, at 
least for x = 1, or that the AF correlations were destroyed 
in the vicinity of the Zn substituents [^4 23 . Also some 
data have been interpreted as due to the total disappear- 
ance of the pseudo-gap in the vicinity of Zn. Finally, the 
dynamics of the local moment |l5| ] appears to be quite 
different in La2Cu04:Al than in our results. 

To clarify the situation, we present in this article an ex- 
tended report of our experimental data, and perform an 
exhaustive comparison with the literature. We examine 
in detail the effect of Zn on 89 Y NMR, in oriented pow- 
ders of YBC0 6+a; :Zn y with 0.5% < y < 4 %, for x = 0.64 
and 1. NMR, being a local probe, provides useful infor- 
mation about the impurity induced short-range and long- 
range effects in the metal via an analysis of the lineshape, 
Knight shift, and linewidth. In section II, we present the 
experimental details regarding sample preparation and 
the procedures adopted for NMR measurements. In sec- 
tion III, a thorough description of the results of our NMR 
work, allows us to highlight the differences in the effect of 
Zn doping in underdoped and overdoped YBCO. In sec- 
tion IV, the NMR shift, linewidth, and relaxation rate 
data are analyzed considering Zn induced local moments 
on the neighbouring Cu. A contrasting comparison with 
the other experimental results introduced in this section 
is contained in the last subsection of IV. In the conclusion 



section we summarize our overall view of the experimen- 
tal situation, and discuss several theoretical works which 
have paid some attention to the induced magnetism in 
cuprates. 



II. EXPERIMENTAL DETAILS 

Samples of YBCOe+x^Zn were prepared by conven- 
tional solid state reaction techniques as described else- 
where [ p4| . Large (single crystal) grain (size > 50 mi- 
crons) samples were made which were finely ground be- 
fore oxygenation (see Ref. |2^| for further details regard- 
ing characterisation). In order to prepare the samples 
with maximum oxygen content, oxygenation was done 
at ~ 300 °C for a long period (> 10 days) which en- 
sured homogeneity of oxygen content. For preparing 
samples with a reduced oxygen content, the maximally 
oxidized samples were treated in vacuum in a thermobal- 
ance at variable temperatures, up to 450 °C. The samples 
were quenched to room temperature when the equilib- 
rium oxygen content was reached. 

In the case of YBCOe+z (without Zn), when the max- 
imum oxidized samples were deoxidized to a point where 
the sample decomposed, the weight loss corresponded to 
5x (= x max - x m j„) = 1.0. However, on addition of Zn, the 
actual maximum value of <5x which could be reached, pro- 
gressively decreases and equals 0.92 for 4% Zn. Since for 
the Zn doped samples, the ortho-tetra structural phase 
transition still takes place for an oxygen content of x m j n 
+ 0.45 (as in YBCOe+z without Zn) (^6), it appears that 
Xmm = 6.0 in the Zn doped samples while x max linearly 
decreases from 7.0 for 0% Zn to 6.92 for 4% Zn. These 
samples with specific oxygen and zinc contents were then 
fixed in Stycast 1266 and cured overnight in a field of 7.5 
Tesla in order to orient the grains with the c axis aligned 
along the applied field direction. 

NMR measurements were performed by standard 
pulsed NMR techniques. We observed the spin echos af- 
ter a 7r/2 — 7r sequence followed by a 37r/2 — it sequence. 
A perfect inversion of the spin-echo in the latter relative 
to the first sequence ensured the correctness of the ir/2 
pulse length (about 13 /isec at room temperature). The 
89 Y shift was measured with respect to a standard YCI3 
solution. 

The 89 Y spin-lattice relaxation time T\ was determined 
using a 7r/2 — 7r sequence, with a repetition time t rep . An 
exponential fit of the nuclear magnetization (obtained 
from a Fourier Transform, FT, of the time domain spin- 
echo signal) as a function of t rep allowed us to deduce T\ 

For YBC06.64 ; Zn, nn resonances are seen (see Fig. 1) 
at low temperatures (T < 150 K). The relative intensity 
of these nn resonances could be enhanced by repeating 
the pulse sequence at a fast rate (t rep ~ 20 sec). Indeed, 
the Ti of the outermost satellite (~ 10 sec), was found 
smaller than that of the main line (~ 100 sec). The fact 
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that the latter has a reduced intensity in such an experi- 
mental condition thus allows us to fix accurately the po- 
sition and the width of the nn resonances. The mainline 
being the narrowest and the most intense, its position 
and width were easily determined with a long repetition 
time, allowing full recovery of the mainline signal. Using 
the positions and widths of the nn resonances determined 
in the manner indicated above, the relative intensities of 
the various lines in the spectrum were determined (by fit- 
ting the lineshape to a sum of three gaussians), for t rep 
> 5Ti of the slowest recovering component, so that all 
the components had fully recovered. The T\ 's of the in- 
dividual lines were determined from an exponential fit of 
their intensities (in the FT) with respect to the repetition 
time. Each T\ measurement took about 15 hours. 

III. RESULTS 

In the following, we present the doping and tem- 
perature dependence of various NMR parameters in 
YBC06+a;:Zn. We shall first report the existence of ad- 
ditional NMR lines detected in the underdoped samples. 
Their characteristics (shift, width, and intensity) enable 
us to associate them with 89 Y nuclei near-neighbours of 
the Zn substituted on the CuC>2 planes (Section III-A- 
1). The results on the main resonance line, which cor- 
responds to 89 Y sites far from the substituted Zn, are 
reported next and compared to those in the pure system 
(Section III-A-2). Spin-lattice relaxation data on the nn 
and main resonance lines are reported in Section III-B. 

A. Resonance line shift and width 

1. Near neighbour resonances 

We will argue here that the additional resonances de- 
tected in YBC06.64:Zn are not seen in YBCOyiZn and 
that the additional resonances are intrinsic and a direct 
effect of Zn substitution. 

As seen in Fig. 1 the nn-resonance positions depend 
on the sample orientation with respect to the applied 
field. Furthermore (see Fig. 2), the relative intensity of 
the outer line increases with Zn content while its posi- 
tion is unchanged. We also see that, YBCCvZn spectra 
(Fig. 3 (a)) do not show additional lines in the tempera- 
ture range of our measurements (80 < T < 350). While 
this is unambiguously evident for the outermost reso- 
nance, the absence of the middle resonance in spectra 
of YBCCviZn is perhaps not immediately obvious. By 
measuring the YBCOyiZn lineshape with a fast repetition 
rate (so that the middle resonance might be enhanced, 
relative to the mainline, due to its shorter Ti), we see 
(Fig. 3 (b)) in fact, that the lineshape of YBCOytZn is 
unaltered by fast repetition (up to one-fourth of T\ of 
the mainline). If there is any change, it is in fact the 



high frequency tail that has a somewhat reduced relative 
intensity, indicating that the tail has a longer T\ . This is 
in keeping with our understanding that the upper tail in 
the lineshape of YBCOy:Zn appears due to those regions 
of the sample which are not fully oxidized and hence have 
a longer T\. In short, YBCC>6.64:Zn has additional 89 Y 
resonances while YBCCVZn does not. In view of the 
abovementioned facts, the additional resonances are not 
due to spurious phases since those should be present inde- 
pendent of the oxygen content (the deoxygenated samples 
are obtained merely by vacuum reduction of YBCOyiZn 
at low T (< 450 °C)). 

In order to identify the origin of these lines we have 
therefore performed quantitative analyses of the spectral 
intensity. Experimental lineshapes for YBC06.64 ; Zn, ob- 
tained with repetition times much longer than the spin- 
lattice relaxation times T\ , were fitted to a sum of three 
gaussians, where the line position and width of the two 
outer lines had been reliably fixed from the short repe- 
tition time spectra. The spectra along with the fits are 
shown in Fig. 4, while the variation of their intensity as 
a function of Zn content is shown in Fig. 5. 

It is of course quite natural to expect that the most 
affected outer line should be associated with the Y nuclei 
nn to the Zn atoms. In the dilute limit, the intensity from 
a purely statistical occupancy of a single neighbouring 
site of Y by Zn, for an in plane concentration c, is 8c(l — 
c) 7 for the 1 st shell (curve A in Fig. 5). But, as Zn 
induces a significant shift of the 2 nd nn Y sites as well, 
we also need to ensure that the 2 nd nn to Y is unoccupied 
by Zn. The corresponding intensity for a purely random 
statistical occupancy would then be modified to 8c(l — 
c) 15 (curve B in Fig. 5) which yields a smaller intensity 
for large Zn concentrations. We see that the intensity 
of the outer line is then consistent with that of the I s * 
nn shell, assuming that all the Zn are substituted in the 
planes ( in Fig. 5 we have taken c = 1.5$/). As for the 
middle resonance, the expression for the intensity due to 
the occupancy of a single Y 2 nd nn site by Zn (with the 
1 st nn unoccupied) is 16c(l — c) 23 in the dilute limit, 
which is much smaller than the experimental intensity. 
If the 2 nd and 3 rd nn are occupied by Zn with the 1 st nn 
unoccupied, the intensity would be (8c(l — c) 7 + 16c(l — 
c) 15 )(l-c) 8 (curve C in Fig. 5). The assignment for the 
middle resonance is not so clear, but for dilute samples 
its intensity is consistent with that of total occupancy of 
the 2 nd and 3 rd nn, with the 1 st nn unoccupied by Zn. 

The T-dependence of the nn line-shifts shown in Fig. 6 
is seen to be Curie-like with a negative hyperfine cou- 
pling. This Curie-like behaviour is usually observed for 
local moments and justifies this denomination that we in- 
troduced in |Q , although the actual magnitude and exact 
origin of this local moment behaviour will only become 
clear hereafter. 

The linewidth of the nn resonances is found to in- 
crease with decreasing T (Fig. 7). The linewidth which 
increases with Zn concentration may be associated with 
the RKKY-like interaction between the Zn induced local 
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moments. With an increase in the concentration of local 
moments, we might expect a frozen magnetic state (most 
probably a disordered spin-glass) at some point in tem- 
perature. An estimate for this is provided by a simple 
analysis in Section IV-E. 

2. Main resonance 

The temperature dependence of the shifts, Aif(T), of 
the mainline for YBC06.64 : Zn is shown in Fig. 8. As re- 
ported before , the mainline shift does not significantly 
depend on the Zn content, and the average carrier den- 
sity at long distance from Zn carrier density is therefore 
nearly unaffected by Zn substitution, at least for dilute 
concentrations of Zn for which the sample is still metal- 
lic. A slight offset with respect to the pure YBC06.64 
is however evident at higher Zn doping levels and might 
be due to a small increase in the carrier concentration. 
Similar slight offsets are detected for the 17 NMR shift 
of these compounds |^7j , but would rather correspond to 
a minute decrease in the hole content. In this latter case, 
a slightly incomplete oxygen loading of the starting sam- 
ples might result from the fact that it has to be achieved 
in sealed vials, and not in a flowing oxygen atmosphere, 
to facilitate 17 enrichment. 

It should be mentioned that in Ref. 0, the main- 
line shift in an unoriented YBCOg.64:Zn4% sample, had 
showed a slight upturn at low temperatures. Since the 
satellite intensities constitute a significant fraction here 
and cannot be clearly distinguished from the mainline, 
the line-position obtained from the peak did not repre- 
sent the true position of the main line. In the present 
work, the different components of the spectra have been 
analysed in the fits with different repetition times so that 
the true position of the main resonance is deduced and 
shows no upturn. 

The width of the mainline for YBC06.64 : Zn l/ (see 
Fig. 9) has a T-dependence similar to that of pure 
YBC06.64 7 in that it initially decreases with decreas- 
ing temperature and then shows an increase below 120 
K which is sample dependent. In the pure system, the 
linewidth can only be associated with a small macro- 
scopic distribution of chain oxygen content (which we es- 
timate of about ±0.02 for most oxygen contents) which 
results in a distribution of shifts at high temperatures. 
The T-dependencies of the shifts around the YBCOe.64 
composition are such that while the shifts are measur- 
ably different at room temperature, their magnitudes be- 
come nearly the same at low-T Q and the NMR line 
becomes narrower. Therefore the width due to a distri- 
bution of oxygen content decreases at low-T. The mag- 
nitude of the width increases with Zn doping, partly due 
to long-distance effects of the spin-polarisation from the 
Zn-induced local moments. The T-dependence of the 17 O 
NMR width is found to be much larger and provides sup- 
plementary information which is analysed in detail by 



Bobroff et al. §7|. 

Turning to the fully oxygenated YBCOyrZn samples, 
we find that here again the mainline shift is nearly inde- 
pendent of Zn concentration (Fig. 10). However, due to 
the broadening of the line at low-temperatures, we cannot 
unambiguously determine whether the maximum seen in 
the T dependence of 89 Y NMR shift of pure YBC0 7 , with 
T(K max ) only slightly larger than T c (and having a pos- 
sible connection to the pseudo-gap) has shifted to lower 
temperatures or altogether disappeared. The Curie-like 
broadening (Fig. 11) in YBCOyiZn is indicative of a dis- 
tribution of magnetic contributions to the line positions. 
This RKKY-like broadening must originate from a mag- 
netic state which develops around the doped Zn. 

For oxygen contents intermediate between O7 and 
6 . 6 , the T-dependence of the 89 Y shift in YBC0 6+2; :Zn 
is qualitatively similar to that in YBCOg+z (Fig. 12(a)). 
However, the sharp decrease in the shift that occurs 
around 100 K for the slightly oxygen depleted samples 
(YBCO6.95 or so) is absent in the Zn doped samples 
where the decrease in the shift is more gradual with T. 
This might again be due to the difficulty in defining ac- 
curately the oxygen content (and therefore the hole con- 
centration for the Zn substituted samples and especially 
for the large 4 % Zn concentration which has been sys- 
tematically investigated. As for the outer resonance, we 
did not perform systematic investigations versus oxygen 
content. It is however clear in the spectra of Fig. 12(b), 
that the low-frequency tail which monitors the position 
of this outer resonance is progressively nearer to the cen- 
tral line when the oxygen content is increased. Further, 
this outer resonance even disappears for x — 0.92 which 
corresponds to the maximum oxygen content for 4 % Zn. 
The NMR shift of the outer resonance with respect to 
the mainline is therefore progressively reduced with in- 
creasing hole content. 

B. Spin-lattice relaxation 

Next, we present the results of spin-lattice relaxation 
measurements for YBCO 6.64 ; Zn. The data were ob- 
tained on the nn lines and the mainline as detailed in 
the previous section. Representative spectra for vari- 
ous values of t rep are shown in Fig. 13. The resulting 
magnetisation recovery for the three lines for the data of 
Fig. 13 are shown in Fig. 14. Exponential fits have been 
found to apply in all cases as illustrated in Fig. 14. 

Taking such data is obviously not straightforward. A 
high enough signal to noise ratio is required, as seen 
for our spectra displayed in Fig. 13. Other publications 
on 89 Y NMR in YBa 2 Cu 4 8 :Zn f|,|§ are completely 
bereft of T\ data, which corroborates the difficulty in ob- 
taining good data. The relaxation rate of the additional 
lines (other than the mainline) is seen (Fig. 15) to be 
strongly enhanced resulting from local moment fluctua- 
tions as is discussed in section IV-A2. The outermost 
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satellite is the most affected which indicates that it must 
result from having Zn as its l st nn. The mainline T\ is 
nearly unaffected which shows that, for dilute concentra- 
tions of Zn, the planar dynamic susceptibility far from 
Zn is unaffected, in accordance with the NMR shift data. 

In YBCO7, Dupree et al. Jl(| measured the effect of Zn 
doping on Ti at room temperature. They found that the 
89 Y spin-lattice relaxation rate was strongly enhanced 
on Zn substitution. However, our data were in complete 
disagreement with theirs. We therefore repeated mea- 
surements on various batches of samples and at various 
temperatures. In all cases we found that the nuclear 
magnetisation recovery fits well to a single exponential 
(Fig. 16) and that the resulting T\ and its T-dependence 
is not significantly different from that of pure YBCO7, 
as can be seen in Fig. 17. We must therefore conclude 
that limited accuracy was responsible for the observation 
done by Dupree et al. jn|. As there are no discernible 
resonances in addition to the main line in these experi- 
ments on YBC07:Zn, the implication is a much weaker 
induced moment in YBCO7, compared to YBCOg.64:Zn, 
in agreement with our bulk susceptibility data |],[l6| • 

IV. ANALYSIS OF THE EXPERIMENTAL 
RESULTS 

A. Local moments in YBC06.64:Zn 

1. nn NMR shifts 

We recall here, that the distinct, well defined reso- 
nances that we have observed in YBCOg.64 correspond 
to Y near neighbour sites of the substituted Zn. The 
Curie-like T-dependence of the position of the first near- 
neighbour line, and the shortening of its Ti at low-T 
are striking experimental evidence of the occurrence of 
Zn induced local moments. The location, spatial ex- 
tent and dynamics of these moments in YBCOg.64 will 
be discussed first. Occurrence of local moments for the 
slightly overdoped composition YBCO7 is also estab- 
lished through the induced long-distance perturbation of 
the host-spin-magnetization. 

The Zn induced local moments are quite clearly located 
in the vicinity of the Zn, and dominantly on the four 
nearest neighbour O or Cu orbitals. In what follows, 
we shall perform extensive comparisons of the 89 Y NMR 
shift with the Zn induced Curie contribution to the spin 
susceptibility (expressed per mole Zn) |||l6) , 

x = 9lL = ^£1 (!) 

Here Na is the avogadro number and p e ff is the ef- 
fective moment. These comparisons will allow us first to 
rule out a localisation of the moments on the O orbitals 
and then to demonstrate that the local moment is dis- 
tributed on the Cu orbitals. Furthermore, assuming that 



the transferred hyperfine couplings are not modified by 
Zn substitution, we will show that our analysis is consis- 
tent with a locally AF state extended over a few lattice 
sites. 

A local moment could also be present around Zn if 
a hole were trapped on the near-neighbour oxygen or- 
bitals. Two quite different physical situations would oc- 
cur, depending on whether the local moment is located 
on the p n or p a orbitals. Since the p„ are directly ad- 
mixed with the Y s orbitals, a strong positive hyperfine 
coupling would result, contrary to our observation of a 
negative Curie contribution to the shift. Therefore, the 
present experiment implies that this shift component can 
only be induced through the oxygen p a orbitals. 

In undoped YBCOg+rr, the Y NMR shift arises from a 
coupling of the Y nuclear spin with the small fraction of 
holes on the 0(2p CT ) orbitals due to their covalency with 
the Cu(3d x 2_ y 2) holes, while the spin-polarization of the 
doped holes themselves is negligible [jjj. This has been 
deduced from the fact that the covalent admixture of the 
0(2p a ) orbital with the Cu(3d x 2_ J ,2) orbital is about 10 
% [ p9[ , which implies that the hyperfine coupling to the 
oxygen holes should be 10 times larger than its coupling 
to the Cu(3dj.2_ a 2) holes. 

Let us first consider the possibility that the Curie con- 
tribution comes from holes localised on the four 0(2p a ) 
orbitals near Zn. The 1 st nn Y site has six O nn which 
we assume are nearly unaffected by Zn and two O mi 
which would exhibit Zn induced Curie magnetism. The 
net 89 Y shift would be written as 

AK? = (6/8) JC + K° + Sf (2) 

where index a refers to a principal direction, Kf is the 
spin shift of the mainline, K" —2C"/T is the Curie con- 
tribution to the spin shift due to its two I s * nn O with a 
moment, and 5™ is the chemical shift. A least- squares fit 
to the data in Fig. 18(a) allows us to extract the two un- 
known parameters Cf and 5f . The chemical shift values 
thus obtained are 5f (6f ) =144 (163) ± 10 ppm. The 
values of Cf are found to be -14000 (-12300) ± 500 ppm 
K for H||c (H||ab). 

Let us compare then these shifts with the macroscopic 
susceptibility Xc with ubK c — 2£/ft/x c /4, if the moment 
is distributed on the four 0(2p <T ) near neighbour orbitals 
to the Zn. Using C M = 9.2 x 10~ 2 emu K/mole Zn [||l|] 
and the Curie term in the shift deduced above, we get 
H/j/ = -1.6 kG. This is of the order of the hyperfine cou- 
pling expected with Cu and nearly 10 times smaller than 
that expected with oxygen. Moreover, we point out that 
a 2 nd nn Y to Zn would not be coupled to the moment on 
the oxygen (a local moment on oxygen is unlikely to be 
spread over more than 4 sites since it would presumably 
arise from hole localisation). Hence, a second line in ad- 
dition to the main line should not be observed, contrary 
to the data from our experiment. One could imagine that 
one has both, localised hole and weakly affected nn Cu. 
But this would require a large susceptibility on the Cu 
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to give the strong shift of the 2 nd nn line compared to 
the mainline. This eliminates the oxygen p a as a possible 
site for local moments. 

If, however, the satellite shift is modelled as coming 
from a hyperfine coupling to the local moments residing 
on the Cu d x 2_ y 2 orbitals, the relevant equation for the 
I s * nn shift in our model is 

A^ = (5/8)JC + iC 



(3) 



^From a fit of the data to this equation, the chemical shift 
values obtained (see Fig. 18(b)), 81 (5f ) =100 (140) ± 10 
ppm, are only slightly different from Sf (<5" b ) =165 (150) 
ppm found in the pure material . The values of Cf are 
found to be -13100 (-11600) ± 500 ppm. This implies a 
hyperfine field Hhf ~ -3.2 kG/Cu which is slightly larger 
than that for the pure material (~ -2 kG). Such a mod- 
ification of Hh f could be attributed to a corresponding 
change of the Cu(3d a; 3_j / 2)-0(2p .) hybridisation due to a 
displacement of the I s * nn oxygen to Zn . Alternatively, 
if the hyperfine coupling stays unchanged, the actual sus- 
ceptibility on the I s * nn Cu, x(l), is larger than % c /4. 
This would imply that further copper ions would have 
a magnetization anti-parallel to the applied field, which 
might be expected if the local moment develops as an 
AF correlated cloud of copper lattice sites. Considering 
this possibility, the Cu 2 nd nn to Zn will bear a small 
negative susceptibility x(2). 

Fig. 19 illustrates schematically the location and the 
orientation of the Zn induced local moments. In such a 
model, the shifts of the I s * and the 2 nd nn Y will be as 
follows; 



AK% = (6/8)K™ + + Kg 



(4) 
(5) 



Here wc can differentiate the hyperfine fields for the 
first and second nn using fisK" — Hhf(l)xi^) an d 
H B K% = H hf (2)x(2). We can then fit AiiTf - (4/8)if s a 
to a Curie term in addition to a constant and likewise for 
Aif 2 Q - (6/8)X s Q . A fit of the observed shifts of the I s * 
and the 2 nd nn Y to these equations yields the following 
values for the corresponding Curie terms; C\ — -14,530 
ppm K, Cf = + 4630 ppm K. The corresponding chemi- 
cal shifts for the I s * and the 2 nd nn Y are found to be 75 
and 144 ppm, respectively. Assuming the same hyperfine 
coupling for I s * nn and 2 nd nn Y, x(2) = — x(l)/3 and 
therefore the macroscopic susceptibility \c — 8x(l)/3. 
Using n B K^ = Hhf(l)x(l), we get a hyperfine field of 
about -2.35 kG which is closer to the value of -2 kG in 
the undoped compound. Although this picture is then 
compatible with the experimental results, the accuracy 
of the data is not sufficient to ascertain its validity. The 
fact that the intensity of the middle resonance cannot be 
assigned solely to the 2 nd nn of Zn implies that further 
near neighbour sites of the Zn should be taken into ac- 
count. More accuracy would also be required to take into 
account the 3 rd nn and try to estimate the size of the AF 
correlated region around the Zn site. 



2. Spin-lattice relaxation 

We next turn to a discussion of the spin-lattice relax- 
ation rate which is expressed as, 



T\T q yH; uj 



(6) 



where A(q) is the coupling to the magnetic fluctuations 
at wave vector q |30| . The O and Y nuclei are at symme- 
try positions with respect to Cu so that A(q^p) = 0, 
and the fluctuations at qAF are filtered at these two 
sites. Consequently, in YBCOg+a;, the T-dependence of 
(TiT)^ 1 for Y and O is different from that of Cu which is 
dominated by the fluctuations at qaf || . On adding Zn, 



the symmetry around the I s 



Y is broken and this Y 



site becomes then sensitive to the magnetic fluctuations 
on the neighbouring copper ions (either intrinsic to the 
pure compound or due to the local moment). Similarly, 
the magnetic fluctuations are no longer symmetric on the 
2 nd nn, so that the enhanced relaxation rate at this site 
is also connected to local moment fluctuations. 

The sharp increase of (TiT) -1 at low-T on the Y nn 
nuclei (much faster than the corresponding variation on 
the 63 Cu nuclei in the pure compounds) is a direct proof 
that the local moment fluctuations are not those of the 
Cu hole spins of the pure compound. The very existence 
of a Curie contribution to the spin susceptibility indeed 
clearly points out that the fluctuation of the Cu hole 
spins in the vicinity of the Zn are much slower than those 
of the pure host. The present data for (TiT) -1 on the 
near-neighbour nuclei are then good proof of the slow 
fluctuations of the local moment. 

In the case of local moments in noble metal hosts, the 
spin- lattice relaxation of host nuclei nearby the local mo- 
ment is totally dominated at low-T by the fluctuations of 
the local moment (the usual Korringa process via conduc- 
tion electrons is somewhat smaller) |38 . The situation 
here is quite similar for the Y I s * nn of the Zn. We can 
therefore consider that this nuclear spin is coupled to the 
2 nn coppers which bear susceptibilities x c (q, w) if we 
neglect the contribution to T\ of the uncompensated Cu 
spin 2 nd nn to Zn (Fig. 19). The relaxation rate at low-T 
is then given by, 



1 



^(^)^ /( l)2E(^^ 



2fc B T / 7r M 2 



(7) 



where jn/je is the ratio of the nuclear and the elec- 
tronic gyromagnetic ratios and Hhf(l) is the hyperfine 
coupling. In the limit w — ► 0, the summation is given by 
Xc(T)t/2tt where Xc(T) is the local moment susceptibil- 
ity and r is the relaxation time of the local moment spin. 
The fluctuation rate of the local moment spin is usually 
made up of two contributions 



1 

Tint 



(8) 
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where the first term corresponds to the single Zn im- 
purity local moment relaxation to the host spin bath, 
for instance through the exchange with the conduction 
electron spins and the second would correspond to fluc- 
tuations due the coupling between the Zn induced local 
moments which depends on Zn concentration. 

For instance, for dilute local moments in noble metal 
hosts 

= &(k B T)(J exP (tF)) 2 (9) 

'ex " 

if a Korringa relation holds (J ex is the coupling of the 
local moments to the band). In that case, the second 
term is t^] = u int /2-K with w? nt = SJf nt zS(S + l)/3h 2 
where z oc c is the number of nearest neigbour spins 
and Ji n t is the conduction electron mediated coupling 
between impurity spins. 

In our case, the T\ values for the near neighbour res- 
onances did not depend markedly on the Zn content, so 
that the single impurity induced relaxation l/r ex domi- 
nates the results. Further, the local moment spin suscep- 
tibility is Curie- like down to low temperatures. A small 
Curie- Weiss correction \ = C/JT + 9) with 9 ~ 4 K is 
observed for YBC06.64 : Zn 4 % |J. This is in agreement 
with the observed spin freezing temperature of about 3 
K in the sample ||. All these results therefore allow 
to conclude consistently that the spin-lattice relaxation 
is dominated by the spin fluctuations of the isolated Zn 
induced local moment. This should then be written as 

1 2k B T flns2TJ fls2 Xc(T)r ex 

J-l h 7e 27T 

Since x c T is constant, 1/Ti is proportional to r ex . In 
the conventional metallic case, where p(c_f) is indepen- 
dent of the temperature, the spin-lattice relaxation rate 
1 /Ti of the host nuclei (near impurities) is observed to 
follow a Curie- like law 1/Ti oc C/T pq|. The present 
case is clearly more complicated since the host metal it- 
self is strongly correlated, which results in the pseudo- 
gap of the static spin susceptibility and in the 1/(T]T) 
behaviour for the 63 Cu. As the local moment posi- 
tions are commensurate with the Cu hole spin system, 
we might anticipate that a similar anomaly might oc- 
cur for l/(r ex T), which scales with T\/T. We have 
therefore plotted in Fig. 20 this quantity for the three 
data points for Y near neighbours to Zn. Although the 
data are clearly insufficient, they suggest a maximum for 
l/(T e xT), quite analogous to that seen for 63 Cu T\. 

B. Comparison with other experiments 

Following our early report ||, some related experi- 
ments have been performed on the impurity substituted 
cuprates. In various cases, the authors have drawn con- 
clusions which are not in complete agreement with our 
results and sometimes disagree totally. These contradic- 
tions are therefore considered in the following. 



1. Gd ESR 

Janossy et al. @ have done Gd 3+ ESR, in 1% Gd sub- 
stituted YBCOe+j, which are in principle quite similar 
to our experiments. Indeed, the Gd electronic moment 
and the 89 Y nuclear spin are coupled to the copper hole 
spins by similar transferred couplings. They find that the 
g-shift of the Gd ESR has the same T dependence as the 
89 Y NMR shift. They have therefore used the Gd probe 
in Yo.ggGd o.oiB a Cu06+a;:Zn samples as well. The Gd 
spectrum should result from a simple scaling with the 89 Y 
NMR spectrum through the ratio of the hyperfine cou- 
plings, as the relative positions of the main and satellite 
lines scale as the ratio of (vHhf) where v is the operating 
frequency and 

{5u/u) = H hfX (11) 

The fact that they did not detect any nn resonances 
might cast some doubt on the meaning of our results. 
However, we insist here that one should also consider 
scaling of the relaxation rates to arrive at reliable con- 
clusions. The relaxation rates scale as 

1/Ti oc (jH hf f (12) 

and therefore contribute quite differently to the broad- 
ening of the spectra. In the case of NMR, the linewidth 
Av is governed by the susceptibility distribution which 
scales as Su, while in ESR, the T\ process is so efficient 
that it contributes significantly to (and even dominates) 
the ESR linewidth. Using our detailed data for 89 Y nn 
NMR, we can easily calculate the expected contributions 
for the Gd nn ESR, through Eq. (11)-(12). Here we 
shall use 0, Gd H hf = 10 Y H hf and j e = 1.6 X 10 4 
89 7. Using the operating frequencies, v — 15.64 MHz for 
89 Y NMR and v = 245 GHz for Gd ESR, we can deduce 
the satellite separation from the main line and the static 
and Ti contributions to the linewidth at 100 K. These 
are reported in Table 1. We find then that the expected 
width for the Gd ESR nn line is at least four times larger 
than its shift, which justifies that the satellite resonances 
are indeed quite difficult if not impossible to observe. 

The other issue is the suggestion by Janossy et al. that 
the susceptibility corresponding to the pure YBCO7 is 
restored at the Cu neighbouring Zn. Since we observe a 
Curie-like increase of the nn line shift in Zn doped sam- 
ples, increasing to values well above the YBCO7 shift, the 
implication is that the hole content near the Zn dopants 
has not been restored to that of undoped YBCO7. 

2. NMR in Al doped La 2 - x Sr x Cu04 

Recently, Ishida et al. jl3) have performed NMR ex- 
periments on La2- a; Sr a ;Cu04 (LASCO), in which non- 
magnetic Al is substituted on the Cu site of the Cu02 
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planes. Although they were unable to detect the nn nu- 
clei of Al, they could directly detect the 27 Al NMR sig- 
nal. They did find that the shift of the 27 Al NMR has 
a Curie component. Since Al itself does not bear a lo- 
cal moment, this observation can only result from a local 
moment which resides either on the nn oxygen or copper 
orbitals which are coupled to the 27 A1 nuclear spin via 
transferred hyperfine couplings. This observation does 
not enable one to decide the location of the local mo- 
ment. However, by analogy with our results, the authors 
have inferred that it is located on the nn Cu orbitals. 
Their data are important as they confirm that a non- 
magnetic substituent induces a local moment in a cuprate 
different from YBCOg+a:- Ishida et al. |r| measured the 
27 A1 NMR shift and T\ which can be compared with the 
corresponding data on the nn 89 Y NMR in YBCOe+a. 
As we stress below, several qualitative differences in the 
experimental results are evident. 

First, they analysed the T-dependence of the shift and 
susceptibility in Al doped LASCO with a Curie- Weiss 
law with a sizeable Weiss temperature (9 m 50 K). It is 
not so clear whether this high value of 9 is also suggested 
by the 27 A1 NMR results since it might be influenced by 
the reference taken for the 27 Al chemical shift. In any 
case, Ishida et al. do not demonstrate whether this large 
9 corresponds to a genuine single impurity effect or if it 
varies with Al content thereby revealing a large coupling 
between the local moments. In underdoped YBCO, we 
never found any indication for such a large deviation from 
the Curie law, neither from NMR nor from susceptibil- 
ity data. A negative value 9 w -30K has been found by 
Monod et al. |31|] in YBCO:Zn by susceptibility measure- 
ments. However recent data on samples with low content 
of impurity phases, by Mendels et al. |J, establish that 
a significant estimate of 9 requires a correct accounting 
of the susceptibility contribution of the pure compound. 
They deduced 9 ~ 4 K for YBC0 6 . 64 :Zn 4% . 

In order to analyse the 27 Al relaxation rate data, Ishida 
et al. use the simple local moment formulation of Eq. 
(7-10) (with different notations). They provided a quan- 
titative analysis of their data in which the value of Ji n t 
as deduced from their value of Ti nt (a microscopic probe) 
is fully consistent with their value of 9 (from bulk suscep- 
tibility, a macroscopic probe). This analysis would then 
appear to be consistent with and support a simple local 
moment picture. We stress here that not only do their 
results differ qualitatively from ours but that an alterna- 
tive interpretation is possible f^j . In their analysis, they 
deduce a relaxation rate 1/r which only varies slightly 
with T . This would be expected if it were dominated by 
the Ti n t term which is expected to be T-independent in 
the classical local moment picture. Our results, on the 
other hand, have exhibited a totally opposite trend with 
the r ex being negligible. We further find that in their T\ 
analysis, Ishida et a/.have taken a local moment Curie 
susceptibility while their own 27 A1 data yielded a large 
Curie- Weiss temperature [9 = 50 K). Introducing the ac- 
tual 1/(T + 50) dependence of x in Eq. 7-10 lead us to 



deduce l/(r ex T) — 2 x 10 10 (sec K) -1 which corresponds 
to a T dependent contribution to 1 /r which is one order 
of magnitude larger than their own result. This would 
yield a rather large value J ex — 0.17 eV which contradicts 
their expectations. 

However, such an analysis yields only a modest mod- 
ification of the T-independent contribution which be- 
comes 1/rint — 5.6 x 10 12 sec -1 , only a factor of 
two smaller than their result which still corresponds 
to a sizeable value of 9, the Weiss temperature. Of 
course, a significant difference between the two systems 
(YBC06.64 an d LASCO) is their doping range. While 
Lai.85Sro.i5Cu04 should be considered close to an opti- 
mally doped material, in YBCOg.64 we are clearly in the 
underdoped regime which usually displays qualitatively 
different properties. Unfortunately we do not have as 
complete results on YBC07:Zn, which would allow for a 
direct comparison between the two systems. 



C. The case of YBC0 7 :Zn 

We have then clearly demonstrated that local moments 
are induced on Zn substitution in YBCOg.64- We have 
not studied samples with oxygen contents other than 
6.64 and 7 in great detail. However, we have seen, from 
measurements on unoriented samples (see Fig. 12(b)), 
that the low-frequency tail of the spectra which is asso- 
ciated with the outer satellite resonance is less shifted 
from the main resonance position for increasing oxy- 
gen content. We can then conclude that the local mo- 
ment value decreases gradually with increasing x. On 
reaching YBCO7, we find that the nn lines have practi- 
cally merged with the main line. This is confirmed from 
magnetisation data on impurity-phase free samples by 
Mendels et al. |) 16 1, which show that the Curie constant 
for YBC07:Zn is about one-sixth that of YBCOe.64:Zn. 
Assuming the same hyperfine couplings, the expected 
first nn line position is shown in Fig. 21. In view of 
the width due to a distribution of oxygen content, it is 
evident that it will be difficult to resolve any extra res- 
onance even for lower Zn contents. Going to lower tem- 
peratures is ruled out as well due to the relatively high 
T c of the samples. 

We did not succeed either in distinguishing the nn res- 
onance from a contrast of relaxation rate with the main- 
line. We shall see here that such a contrast is not ex- 
pected if we estimate the relaxation rate for the outer- 
most resonance by scaling the YBC06.64 ; Zn data at 100 
K. The Curie term in YBCOyiZn is about one-sixth that 
in YBC06.64:Zn and the density of states at the Fermi 
level /o(cf) can be estimated about 3 times higher, from 
the 89 Y NMR shift data. Assuming that the local mo- 
ment to band coupling J ex stays unchanged, Equations 
9 and 10 allow to deduce a contribution of the local mo- 
ment fluctuations to 1/Ti of « 0.0018 sec -1 , which is 
much smaller than the observed rate in undoped YBCO7 
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(0.03 sec" 1 ). This confirms that for YBC0 7 :Zn, local 
moment fluctuations are indeed difficult to detect on the 
nn Y site. 

In fact, the occurrence of a local moment in YBCOyiZn 
was evidenced first JtJ from the presence of the oscillating 
long distance RKKY spin polarisation of the host Cu 
spins. This was established from the Curie like increase 
of the 89 Y NMR linewidth observed in YBC0 7 :Zn, and 
is confirmed in the present experiments as well as from 
17 NMR linewidth M data. 



D. About AF correlations near the Zn impurities 

The experimental results on Gd ESR in YBCO and Al 
NMR in Al doped LASCO have been interpreted along 
quite different lines. For Janossy et al., the absence of a 
Curie term for the ESR line shift led them to conclude 
that there was no local moment. Instead, they considered 
that the susceptibility of the YBCO system is restored 
near the Zn, as they found an increase of the ESR shift 
with decreasing T. On purely experimental grounds it 
is not clear whether the detected ESR signal involves all 
the Gd spins. We have seen above that the outermost nn 
resonances are not expected to be resolved in the ESR 
data. However, the inner resonance might contribute to a 
wing in the signal, and might explain the observed shift. 

In any case, the present detailed 89 Y NMR data 
demonstrate that the central line is not shifted at all, 
so that the susceptibility is unmodified at a few lattice 
distances from the Zn impurity. Second, the Curie-like 
increase of the shift for the nn reaches values well above 
those observed for pure YBCO7, which implies that the 
hole content near the Zn dopants has not been restored to 
that of undoped YBCO7. Finally, the bulk susceptibility 
data (measured using a commercial SQUID magnetome- 
ter) of Mendels et al. || indicate that the local moment 
susceptibility increases down to 10 K, so that there is 
no doubt about the occurrence of a Curie contribution. 
The fact that the nn lines broaden stongly with decreas- 
ing temperature is sufficient to explain that the Gd ESR 
picks up only a part of the Gd signal which saturates at 
low- T . This is somewhat reminiscent of the situation 
which prevailed in the preliminary 89 Y NMR measure- 
ments done for large Zn concentrations in YBCO [Q . In 
that case, the nn resonances were not resolved and an 
apparent shift of the 89 Y NMR signal was observed. As 
the broadening of the nn lines is expected to be larger in 
the Gd ESR, the measured shift involves a contribution 
of those nn sites and is much smaller than that expected 
for the I s * nn signal. 

As for the analysis of the Al NMR data in LASCO, the 
authors of course do not question the existence of a local 
moment. But, they still consider that AF correlations 
are reduced near the impurity, and that the induced mo- 
ments on the four copper near neighbours are decoupled. 
Further, they even anticipate that a state nearer to that 



observed in the overdoped material prevails at distances 



just greater than the I s 



distance ]13[. However, we 



feel that there is no experimental evidence for such a pos- 
sibility in their work on the LASCO system. The only 
argument advanced by Ishida et al. to support this hy- 
pothesis is the independent experimental evidence found 
in their group for two T\ components in their 63 Cu NQR 
measurements in Zn doped YBCO7, and in YBa2Cu40g 
p3|p^| . In both cases they find that the long T\ compo- 
nent is longer than that observed in the pure system, and 
they therefore associate it with Cu nuclei near the Zn im- 
purities. This leads them to conclude that the magnetic 
fluctuations near the Zn impurites have been suppressed. 
However, the relative magnitude of the two components 
in terms of number of sites and their dependence on Zn 
doping which could support this interpretation has not 
been studied in detail. Most importantly, the underlying 
idea seems to us to contain an essential contradiction. 
Indeed if the AF fluctuations around the Zn were sup- 
pressed, this would imply that Zn is in a classical metal- 
lic environment, which would be totally inconsistent with 
the occurrence of a local moment |32| . 

In the superconducting state, Ishida et al. do find 63 Cu 
NQR relaxation rates much larger than those found in 
the pure system, which indicates the existence of states 
in the superconducting gap. This is also seen from Yb 
Mossbauer experiments on samples in which a small frac- 
tion of Y has been substituted by Yb 35 1. States in the 
gap in the superconducting state induced by Zn impu- 
rities were also seen by neutron scattering experiments 
p9[ . Those states are found at a scattering vector of (tt, 
7r), even for YBCOyiZn , while a scattering at (it, it) for 
this oxygen content can hardly be detected in the pure 
system. These experiments are thus direct evidences in 
favor of the persistence of AF correlations in the vicinity 
of the impurities. 

All these observations support the main point which 
we have been advocating, that the AF correlations are 
at least maintained, and perhaps even strengthened near 
the Zn impurities. In such a case the local moment can- 
not be considered as formed of four independently fluctu- 
ating moments on the four Cu sites nn to Zn, but rather 
as an extended state involving further neighbours, and in 
which the Cu nn to Zn are ferromagnetically correlated 
and fluctuate as a single entity. 

We therefore think that the experimental observation 
done by Ishida et al. on the normal state 63 Cu NQR T\ 
might have a quite distinct interpretation. A more sys- 
tematic study, possibly with different impurities might be 
needed to clarify the origin of the longer T\ component. 

In conclusion, it seems to us that the existing exper- 
iments do not contradict the main point of view origi- 
nally proposed, i.e. that the local moments induced by 
Zn are associated with the correlated nature of the Cu02 
planes and that AF correlations might even actually be 
enhanced around Zn. 







E. Induced spin polarisation at large distance from 
the Zn 

Up to now we have mainly considered the magnetic 
moments induced near the Zn impurities. In noble met- 
als hosts, any local charge perturbation is known to in- 
duce long distance charge density oscillations (also called 
Friedel oscillations). Similarly a local moment induces a 
long distance oscillatory spin polarisation (RKKY) which 
has an amplitude which scales with the coupling J ex of 
the local moment with the conduction electrons. This 
oscillatory spin polarisation gives a contribution to the 
NMR shift of the nuclei which decreases with increasing 
distance from the impurity. In very dilute samples, if 
the experimental sensitivity is sufficient, the resonances 
of the different shells of neighbours to the impurity can 
be resolved |38[ . These resonances merge together if the 
impurity concentration is too large, which results then in 
a net broadening of the host nuclear resonance. 

Here, the occurrence of the local moment, even induced 
by the non-magnetic substitution is also a local magnetic 
perturbation in the correlated host. One therefore ex- 
pects a response which will extend to long distances from 
the impurity. Such contributions to the NMR linewidths 
have been found in our work. We shall consider here in 
turn the case of YBCOg.6 an d that of YBCO7. 

1. YBC0 6 . 6 

Indeed, both the central 89 Y line as well as the near 
neighbour resonances have been found to be broadened 
in YBCC>6.6- As seen in figures 7 and 9, these linewidths 
increase at low-T and also increase with increasing im- 
purity content. The central line broadening is unfor- 
tunately only a small fraction of the pure compound 
linewidth, and the temperature dependence of the impu- 
rity induced contribution cannot be extracted accurately. 
Experiments have therefore been performed by Bobroff et 
al. [^7j on 17 nuclei in substituted samples. The larger 
hyperfine coupling of the 17 nuclei with the planar Cu 
as compared to that of 89 Y lead consequently to broad- 
enings of the 17 NMR width, which have been studied 
in great detail both for Ni and Zn substitutions. It has 
been found that in both cases the broadening increases 
much faster than 1/T at low temperatures, contrary to 
what one might expect in a non-correlated metallic host. 
This fast increase is a signature of the anomalous mag- 
netic response of the host which displays a peak near the 
AF wavevector (n, 7r). 

In the present experiments, the broadenings of the I s * 
nn line (Fig. 7) are somewhat related to this long dis- 
tance polarisation induced by the Zn impurity. The large 
increase of the nn linewidth with increasing Zn concen- 
tration is due to the distribution of susceptibility of the 
moments associated with their mutual interaction. In a 
molecular field approach, the Curie contribution K c to 



the shift of a 89 Y nn of a given Zn atom is proportional 
to x(i?o + H m ), where x is the single impurity dimcn- 
sionless susceptibility (— Ci mp /T) and H m the molecular 
field at the moment site induced by other Zn moments. 
This molecular field scales with the magnetization of the 
local moments (H m — kM) and therefore varies as 1/T. 
The linewidth is then related to the root mean square 
value of the molecular field 8H m . Consequently, the 
linewidth due to interaction between the local moments 
(which scales with x$H m ) should scale as 1/T 2 . We have 
therefore plotted in Fig. 22 the quantity T 2 AH corr /H 
versus T, where AH corr = AH nn — AH pure is the in- 
crease of the nn satellite linewidth (Full Width at Half 
Maximum) with respect to that of 89 Y in pure YBCC>6.6- 
We can see that T 2 AH corr /H is nearly T-independant as 
expected from such a simple model. Let us point out that 
H m should in principle behave as the long distance spin 
polarisation detected by 17 O NMR, and should then in- 
crease faster than 1/T at low temperature. Although the 
experimental accuracy on the NMR width is not great, a 
large increase of T 2 AH corr /H is not observed at low T. 
More detailed and possibly more accurate experiments 
are required to better understand whether other contri- 
butions to the nn linewidth have to be considered as well. 

^From our data we can however get an overestimate 
of SH m from a comparison of the magnitude of the 
linewidth with the actual shift of the nn line. Assuming a 
gaussian shape for this resonance AH corr /2.36 is simply 
proportional to x8H m , while the shift K c Hq is propor- 
tional to xHq. Therefore, 5H m — 0.42 AH corr /K c . Fur- 
ther, from the analysis of Eq. 3, K c ~ 0.024/T. From the 
discussion above, AH corr = 2.36kc 2 mp H /T 2 and from 
Fig. 22 for 1 % Zn, AH corr ~ 2H /T 2 . Then, for an 
applied field Hq — 7 Tesla, we deduce 

SH m = 250/T(Tesla/%Zn) (13) 

The molecular field becomes comparable with the ther- 
mal energy for ksT = /j, e ffH m , which for a measured 
Hefj ~ 0.8/i_b in YBCO6.6 corresponds to about 1.2 
Tesla/K. Therefore the temperature at which a spin-glass 
freezing of this system should occur can be estimated 
to be about 15 K for 1 % Zn. This number deduced 
from this rough analysis is somewhat higher than that 
obtained from the Weiss temperature measured by static 
susceptibility by Mendels et al. M, which does not ex- 
ceed 4 K for 4 % Zn. Apart from the above mentioned 
possible experimental limitations, this difference could be 
linked with the fact that we are dealing here with a 2D 
Heisenberg spin system, for which quantum fluctuations 
reduce the spin-glass ordering temperature to T g = 
p6| . A finite value for T g would then only result from 
weak interplane exchange couplings. 



10 



2. YBCOi 

The broadening has been found to increase as 1/T at 
low-T, for this slightly overdoped composition for which 
the planar susceptibility of the pure system has little In- 
dependence. This increased linewidth at low-T is a di- 
rect proof of the existence of a local moment behaviour 
induced by Zn for this overdoped system fl38[ . We have 
seen that the absence of near neighbour resonance lines 
is also an indication that the effective moment is very 
small, which confirms the susceptibility data of Mendels 
et al.. A similar observation has been made by Bobroff et 
al. [§7| through 17 NMR data. Initially such a broad- 
ening can be explained with the RKKY-like broadening 
induced by the local moments. However, in the under- 
doped case it has been shown that the response of the 
correlated electronic system is quite distinct from that 
of a free electron g the 17 NMR linewidth ex- 

hibits an anomalous T dependence The NMR data for 
the YBCO7 composition, both for 17 and 89 Y, does 
not display such an anomalous T dependence, and one 
might wonder whether a conventional RKKY broadening 
is then recovered. Such an approach has been used in our 
initial report @. However, the very fact that a T- depen- 
dent magnetic behaviour is induced by Zn substitution 
is an indication that the correlated nature of the elec- 
tronic state has not disappeared in YBCO7. This is also 
established by the well known anomalous non-Korringa 
T-dependence of 63 Cu p^ . Therefore a direct test of the 
shape of the spatial dependence of the impurity induced 
spin polarisation should give us information on the im- 
portance of these correlations. This aim will be pursued 
in the future with careful studies of the NMR lineshapes, 
which are expected to be more sensitive to the detailed 
shape of the spin polarisation [p7| . Such experimental 
studies will be undertaken on 17 NMR which possesses 
a larger signal to noise ratio. A comparative discussion 
of the induced spin polarisation as sensed by the 89 Y and 
the 17 nuclei will therefore be performed in the future. 



V. CONCLUSIONS 

A large variety of conclusions have been drawn and 
various questions have been raised from the present re- 
sults. They address different points extending from the 
materials properties to detailed questions on the elec- 
tronic structure of the impurities and their influence on 
superconductivity. 

First, concerning the physical chemistry of the 
cuprates, the intensity of the near neighbour resonances 
allowed us to calibrate the amount of Zn substituted on 
the planar Cu sites. Our result is the strongest experi- 
mental proofthat the Zn substitutes dominantly on this 
planar site, up to 3 % Zn, and within 10 % experimental 
accuracy. 



We have confirmed the influence of Zn impurities on 
the phase diagram of the cuprates in the underdoped 
regime. The implication that the static and dynamic 
susceptibility far from the impurity is unaffected by Zn is 
borne out by our shift and relaxation data. This demon- 
strates that the related q = pseudo-gap is not modi- 
fied. The change of the macroscopic susceptibility is only 
associated with modifications of magnetic properties in 
the vicinity of the impurity. Kakurai et al. |lq ] initially 
suggested, on the basis of their neutron scattering exper- 
iments, that the pseudo-gap vanishes at q =(tt,tt) while 
the gap at other q values is unchanged. However, the 
neutron data of Sidis et al. (T?^] in fact suggests that the 
pseudo-gap at q =(w, n) does not vanish but that some 
states appear in the pseudo-gap. Those could also be as- 
sociated with the local magnetic modifications induced 
around the Zn. In a scenario in which the pseudo gaps 
would be associated with the formation of local pairs at 
high-T, these results ndicate that impurities do not pre- 
vent the formation of local pairs except possibly in their 
vicinity. 

What are the actual magnetic properties in the vicin- 
ity of the Zn impurity? Although our early experiments 
had given strong proofs of the occurrence of a local mo- 
ment behaviour induced by non-magnetic Zn impurities, 
the validity of this observation has been periodically put 
into question. The significance of the nn 89 Y NMR re- 
sults has been, for instance, questioned from the absence 
of detectable nn resonances of Zn in the ESR experi- 
ments on Gd/Y substituted underdoped samples. We 
have clearly shown here that the large expected relax- 
ation rate induces a broadening of the Gd ESR nn lines 
which prohibits their detection. The authors have also 
concluded from those Gd ESR experiments that the full 
density of states corresponding to pure YBCO7 is re- 
stored near the Zn impurity. The fact that the 89 Y I s * nn 
resonance is found to display an NMR shift much larger 
than that of the optimally doped compound at low- T, is 
clear evidence against this idea. On the contrary, the sus- 
ceptibility of the Cu nn to Zn is found to present a Curie 
like T-dependence, hence the "local moment" denomina- 
tion, which we have been using throughout. This local 
moment behaviour is confirmed in YBCOg.6 by macro- 
scopic susceptibility SQUID data |) 16 1. 

It is clear that the observed local moment behaviour 
is original inasmuch as it is the magnetic response of the 
correlated electron system to the presence of a spinless 
site. The perturbation induced by Zn extends at least 
to the four nn copper sites, but we have shown that, 
in underdoped YBCOg.6, our data are compatible with 
a local dynamic AF state which extends over more Cu 
sites. Although the present NMR data are not sufficient 
to allow us to determine the actual extension of this state, 
the width of the neutron scattering peak at (tt, it) which 
is found to develop at low-T within the pseudo-gap in 
presence of Zn |l^] , corresponds to a real space extension 
of at least 7 A. 

Various theoretical arguments in favor of the occurence 
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of a local moment in presence of a spinless site in a cor- 
related electronic system have been advanced |],^9] ElJ . 
As complete understanding of the magnetic properties of 
pure cuprates is far from being achieved, it is no surprise 
that the present theoretical descriptions of the impurity 
induced magnetism are rather crude, and for example, 
do not address its microscopic extent. Our results might, 
however, be put in parallel with recent theoretical work 
on undoped quantum spin systems. For instance Mar- 
tins j|D predicts static local moments induced due to 
doping 5=1/2 Heisenberg AF chains or ladders with 
non-magnetic impurities. NMR experiments on the S 
= 1/2 Heisenberg chain system Sr2Cu03 are consistent 
with the prediction of an induced local moment with a 
large spatial extent along the chain |^| . In this undoped 
insulating quantum liquid, the response is then purely 
magnetic. Since the parent compound to YBCO super- 
conductors is a 2D Heisenberg AF and dynamic AF cor- 
relations appear to persist even in the metallic compo- 
sitions, appearance of local moments on many Cu sites 
near to the doped Zn might well be anticipated. 

In the slightly overdoped YBCO7, the local moment 
could initially only be detected through the induced long 
distance spin polarisation H]. A local moment induced 
by non-magnetic Al substituted on Cu is also detected 
in optimally doped LASCO from 27 Al NMR . The fact 
that we could not resolve the nn signal in YBCO7 is 
consistent with the weak magnitude found for the Curie 
like contribution to the local susceptibility. The decreas- 
ing magnitude of the moment with increasing hole doping 
could carefully be monitored by direct SQUID measure- 
ments [lq| . This decrease could be linked experimen- 
tally with a decreasing screening radiusby the conduction 
band. However, the magnetic states which are detected 
within the spin-gap at low-T by neutron scattering JToj ] 
exhibit a short magnetic correlation length, so that the 
spatial extent of the local moment also decreases with in- 
creasing hole doping as does the AF correlation length in 
the pure system. Altogether, our experiments cannot, at 
present, distinguish the respective roles of the screening 
radius and the AF correlation length in defining the local 
moment magnitude and spatial extent. 

Another important question which arises then concerns 
the coupling of the defect local moment to the host. For 
magnetic impurities in simple metals, an exchange cou- 
pling J ex between the local moment and the conduction 
electron spins usually occurs, and determines some of the 
thermodynamic properties of the local moment. For in- 
stance the fluctuation rate of the local moment (1/r) is 
directly determined by J ex , and follows a Korringa re- 
lation in classical metals. This fluctuation time can be 
estimated from nuclear spin lattice relaxation data. In 
the YBCO system we could only obtain such measure- 
ments in the underdoped regime on the 89 Y nn of Zn. 
From these we could show that only weak contributions 
to l/Xi are expected on the 89 Y nn in the optimally 
doped case, and could not be sensed within experimental 
accuracy. Direct measurements of the 27 A1 T\ are on the 



contrary sensitive enough in the optimally doped case in 
LASCO, as seen by Ishida et al. Their results, although 
they establish the occurrence of a local moment induced 
by the spin- less Al 3+ substituent, differ markedly from 
those obtained in YBCO. A large temperature indepen- 
dent contribution to the shift and local moment fluctua- 
tion time is detected, contrary to our observations. The 
origin of the difference might be 

i) linked with the larger valence of Al 3+ , i.e. the charge 
difference with respect to host planar Cu 2+ , 

ii) a peculiarity of the LASCO system, as indeed the 
physical properties of this system do not appear to fit in 
a universal picture with the other cuprates (see Ref . Q ) , 

iii) or merely a difference between underdoped and op- 
timally doped systems, as the experiments could not be 
performed on the two systems under similar conditions. 

Further experiments will certainly permit to make a 
decision between these possibilities. Currently, experi- 
ments do not permit a clear indication on the applica- 
bility of an exchange model. In conventional metallic 
systems, the local moment couples through J ex to the 
electron bath and an oscillatory RKKY polarization oc- 
curs in the band. Therefore J ex can be usually estimated 
from the broadening of the host NMR |2^Q. Applying 
the standard RKKY theory yields values of the exchange 
coupling which are very large [pi. But, we have recently 
shown in Orsay [^7) that, at least in the underdoped 
regime, the behaviour of the 17 linewidth does not fol- 
low the expected RKKY T dependence at all, i.e. the 
NMR width does not scale solely with the impurity mag- 
netization. Let us note here that whatever the method 
used ||, (l3), the estimates of the coupling constant 
are presently such that if one applies a simple exchange 
model, one would expect a large Kondo temperature Tr 
and correspondingly, a spin susceptibility which would 
deviate from the Curie dependence at T ~ Tr and satu- 
rate below. From SQUID data, Mendels et al. concluded 
that this is not the case, and that Tk does not exceed 
at most a few K, in the underdoped YBCO compounds. 
Such a Kondo-like effect was a candidate mechanism for 
the reduction of the magnitude of the local moment in 
YBCOyiZn (sec for instance Ref. [Q). But obviously 
the Kondo model needs to be revised in the context of a 
strongly correlated electron system. Such difficulties had 
been already pointed out by Hirschfeld Jl7| in view of our 
preliminary experimental results. 

In conclusion, we have detailed here the experimental 
evidence for the occurence of a local moment behaviour 
induced by spinless substitutions on the Cu site in Cu02 
planes of cuprates. The existence of original magnetic 
behaviour induced by non-magnetic substitutions can be 
anticipated from current theoretical treatments of un- 
doped low-dimensional spin systems. However, the de- 
tailed experimental observations reported here on doped 
cuprates do not have a thorough interpretation from the 
theoretical standpoint. We suggest that further experi- 
mental and theoretical efforts regarding these properties 
are essential to lead us towards a comprehensive descrip- 
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tion of the magnetic and superconducting properties of 
the cuprates. 
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Figure Captions 

FIG. 1 89 Y NMR lineshape at 130 
K in ¥6006.64:2)1!% when the duty cycle of the pulse 
sequence, t rep , is 20 sec, with the sample c axis aligned 
parallel or perpendicular to the applied field. The rela- 
tive intensity of the nn lines is enhanced here, since they 
have a T\ comparable to t rep while the mainline T\ is 
much longer. 

FIG. 2 

89 Y NMR lineshapes at 100 K in YBCOe.^Zn^, for 
the sample c axis aligned parallel to the applied field H. 
The relative intensities of the outer and middle lines are 
seen to qualitatively increase with increasing y. 

FIG. 3 (a) 89 Y NMR lineshape at 90 K in YBC0 7 :Zn 1% 
showing absence of resolved nn lines indicating that the 
induced local moment magnitude is weak in YBCOviZn 
as compared to that in YBCC>6.64:Zn. Also shown is the 
decomposition of the lineshape for YBC06.64:Zn 1 % into 
three gaussians. 

FIG. 3 (b) 89 Y NMR lineshape at 80 K in 
YBC07:Zn 2 % for different repetition rates. The arrow 
indicates the position of the middle line in YBC06.64 : Zn. 
The lineshape stays nearly unchanged indicating the ab- 
sence of any components relaxing faster than the main- 
line. 

FIG. 4 Fully relaxed 89 Y spectra in YBC0 6 . 6 4:Zn y% . 
The solid lines are fits to three gaussians as explained in 
the text. (a)l% Zn (b)2% Zn (c)4% Zn. 

FIG. 5 Variation of the fractional nn line inten- 
sity (integrated) as a function of Zn content y for 
YBCOg^^Zn^. The solid lines correspond to variations 
as expected from statistical models as explained in text. 
The intensity of the outermost line is seen to be in near 
agreement with that expected from a 89 Y nuclei nearest 
to the doped Zn. The middle line intensity might agree 
with that expected from a combination of 2 nd and 3 rd 
nn 89 Y nuclei. 

FIG. 6 89 Y nn line shifts K versus temperature T for 
YBC06.64 : Znj,%. The outer line shift is seen to have a 
strong upturn with decreasing T indicating a coupling to 
a Curie-like susceptibility. The solid lines arc drawn as 
guides to the eye. 

FIG. 7 Linewidths normalized to the applied magnetic 
field AH/H of the nn lines for YBCOg^Zn^ are seen 
to increase with decreasing T. The solid lines are drawn 
as guides to the eye. 

FIG. 8 89 Y mainline shift K versus temperature T 
for YBC06.64:Zn y % is nearly unchanged from that of 
YBC06.64 indicating little change in the hole content 
with Zn doping. 

FIG. 9 Normalized linewidth AH/H of the mainline 
versus the temperature T for YBC06.64 : Zn y %. Unlike 
YBCO7 which has a nearly T independent linewidth, 
YBC06.64 linewidth decreases below 120 K and increases 
again at much lower temperatures. Qualitatively, the 
linewidths for Zn doped YBCOg.64 ar< 3 higher than the 
undoped compound. The solid lines are drawn as guides 
to the eye. 
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FIG. 10 The T variation of 89 Y mainline shift K for 
YBC07:Znj,% is unchanged from that in YBCO7 again 
indicating that the hole content is nearly unchanged with 
Zn doping. Note that the chemical shift reference is 
about 150 ppm, hence the change in susceptibility on 
Zn addition is at most 4 % of the susceptibility. 

FIG. 11 Normalized linewidth AH/H of the main- 
line versus the temperature T for YBCC>7:Zn y %. The 
linewidth increases in a Curie-like manner with decreas- 
ing T. Also, the linewidth is larger for larger Zn contents. 
Solid lines are drawn as guides to the eye. 

FIG. 12(a) 89 Y mainline shift K versus temperature 
T for YBCO6+x:Zn4o/ . The T-dcpcndcnce is similar to 
that of YBC0 6+a; . 

FIG. 12(b) Variation of lincshape as a function of x 
for YBCC>6+iE:Zn4%. The low frequency tails which cor- 
respond to the shifted satellite lines are seen to appear 
already for oxygen content x = 0.84. 

FIG. 13 Spectra for YBC0 6 .64:Zn 1 % for various t rep 
values. It is clear that the outer line recovers its full 
intensity for much smaller t rep values than the mainline 
and hence has a much shorter T\ than the mainline. 

FIG. 14 Analysis of the relaxation rate data of Fig. 13 
for YBC06.64:Zn 1 %. The data have been fitted to a 
sum of three gaussians and the deduced intensities corre- 
sponding to the main and the nn lines are plotted versus 
the repetition time of the pulse sequence t rep on a semi- 
log scale. The solid lines are fits to a single exponential 
recovery. 

FIG. 15 (a) 89 Y nuclear spin-lattice relaxation rate di- 
vided by temperature \jT\T versus temperature T for 
the mainline and the middle satellite in YBC06.64 : Zn 2/ %. 
(b) these data are for the outermost satellite line. The 
nn lines are seen to have a shorter T\ than the mainline. 

FIG. 16 Nuclear magnetization corresponding to the 
mainline versus the repetition time of the pulse sequence 
t rep for YBCOy:Zn y o/ . The fact that the data can be fit 
to a single exponential (solid line) indicates the absence 
of any other components to the relaxation. 

FIG. 17 89 Y nuclear spin-lattice relaxation rate di- 
vided by temperature 1/T\T versus temperature T for 
YBC07:Zn 2/ %. Note the magnified scale of the y axis. 
As for YBC06.64:Zn, the T\ of the mainline is not af- 
fected. 

FIG. 18(a) Curie component (in addition to a con- 
stant) of the 89 Y shift K versus temperature T for the 
outermost line in YBC06.64 : Zn!%. The solid line is a fit 
(see Eq. (2)) assuming that the local moment is on the 
oxygen atoms, (b) The solid line is a fit (see Eq. (3)) 
assuming moments on the copper atoms. 

FIG. 19 A schematic of the location and the orientation 
of the local magnetization around a Zn impurity. 

FIG. 20 The ratio of the relaxation rate of the local 
moment spin 1/r to temperature T for the I s * nn 89 Y 
nuclei in YBC06 64 ; Zn. r has been determined from Eq. 
(10). 

FIG. 21 89 Y NMR spectrum for YBC0 7 :Zn y% . Ex- 
pected position of outermost line, based on the macro- 



scopic susceptibility data (see text) is indicated by an 
arrow. The actual linewidth even in the pure compound 
is such that any feature at this position cannot be re- 
solved. 

FIG. 22 Variation of the product of the square of tem- 
perature T 2 and the nn linewidth AH corr (linewidth 
of the pure compound has been subtracted from 
the measured linewidth) with the temperature T for 
YBC0 6 . 64 :Zn. 

TABLE I. 89 Y NMR measured values (at about 100 K) 
of the separation of the l si nn line from the mainline 8v, 
the spin- lattice relaxation rate of the 1 st nn 1/Ti, and the 
linewidth of the l si nn Av are listed along with calculated 
values (see text) of the same for Gd ESR in YBC06.64:Zn. 
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